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摘要 : ”近年 来 ， 随 着 GNSS (global navigation satellite system) 技 术 的 发 展 ， 利 用 地 基 和 空 基 
GNSS 电离 层 探 测 的 空间 特征 ， 形 成 了 更 加 完善 的 三 维 电 离 层 观测 系统 ， 获 得 了 丰富 的 多 源 
GNSS 电离 层 观 测 数据 ， 促 进 了 电离 层 在 理论 模型 、 经 验 模型 及 同化 模型 方面 的 进展 ， 同 时 ， 通 
过 结合 中 性 大 气 及 等 离子 层 等 方面 的 研究 ， 为 极端 气象 、 太 阳 活 动 及 磁 暴 等 事件 的 监测 与 物理 
> 过 机 制 研 究 提供 了 数据 基础 与 研究 手段 。 介 绍 了 多 模 地 基 及 星 载 GNSS 电离 层 观测 数据 发 展现 状 ， 
es 并 讨论 了 空 基 GNSS 电离 层 观 测 数据 在 电离 层 顶 部 区 域 经 验 模型 改进 方面 的 进展 。 指 出 电离 层 
^ 经 验 模型 可 作为 同化 模型 的 初始 场 ， 而 电离 层 理论 模型 可 将 其 作为 卡尔 曼 滤 波 的 动力 学 部 分 ， 以 
通过 同化 算法 利用 多 源 数据 重建 或 预测 电离 层 ， 并 阐述 了 多 源 及 多 模 GNSS 电离 层 观测 数据 的 
增加 背景 下 电离 层 同化 模型 构建 过 程 中 新 的 要 求 与 发 展 。 最 后 ， 提 出 自 适应 格 网 与 机 器 学 习 可 成 
CN 为 电离 层 同 化 算法 方面 的 新 研究 方向 ， 以 促进 三 维 电 离 层 研究 。 

Qo x $2 词 ，GNSS， 电 离 层 ， 经 验 模型 ， 同 化 模型 
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作为 日 地 空间 环境 重要 组 成 部 分 ， 电 离 层 对 现代 无 线 电工 程 系 统 和 人 类 的 空间 活动 ( 例 
如 电离 层 与 卫星 导航 系统 定位 的 误差 及 地 震 异常 密切 相关 "”) 有 着 重要 的 影响 。 研 究 电离 
层 不 仅 有 利于 认识 电离 层 本 身 、 寻 找 克 服 电离 层 可 能 造成 灾害 的 途径 和 探求 利用 电离 层 为 
人 类 造福 的 方法 ， 也 有 助 于 推动 地 球 科学 领域 相关 电离 层 理论 和 应 用 问题 的 研究 与 发 展 。 近 
年 来 ， 随 着 全 球 定 位 系统 (global positioning system, GPS) 和 各 国 导 航 星座 的 建成 及 区 域 变 
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] 于 电离 层 研 究 的 新 时 代 。 


GNSS 卫星 星座 的 增加 与 完善 ， 全 球 及 
收 机 的 广泛 使 用 ， 使 得 地 基 GNSS 电离 层 观测 数据 数目 大 量 增加 。 
数据 分 布 与 地 基 监 测 站 分 布 密切 相关 ， 在 海洋 区 域 、 南 美洲 、 
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则 环境 监测 项 目的 开展 ， 迎 来 了 地 基 GNSS 及 空 基 GNSS 等 多 源 电离 层 观 测 数据 应 


局 部 区 域 监测 网 的 构建 ， 以 及 支持 多 模 GNSS 接 
地 基 GNSS 电离 层 观 测 


FERRE R 


， 由 于 地 | 


站 的 


稀 朴 导致 这 些 区 域 的 电离 层 精度 较 差 ， 导 航 


15 


运动 极 慢 ， 即 使 提高 数据 采样 率 ， 对 穿刺 信 


恩 几 何 分 布 的 影响 极 小 ， 数 据 空间 分 辨 率 的 提高 也 并 不 显著 。 由 于 低 轨 卫星 运动 速度 比 导 航 


卫星 快 ， 它 可 以 高 动态 地 进行 全 球 范 
取 的 斜 向 路 径 电 子 总 量 (total electro 
建 。 增 加 低 轨 卫星 电离 层 观 测 数据 ， 可 以 避免 海洋 数据 缺失 。 现 今 


的 低 轨 卫星 3 
的 电离 层 全 球 观测 系统 。 


围 的 电离 层 探测 。 低 轨 卫 星 一 般 可 以 获得 其 定位 数据 提 
n content, TEC) 及 电离 层 掩 星 观测 用 于 电离 层 模型 构 


由 于 电离 层 是 
时 间 积 累 而 误差 增 大 ， 而 


误差 。 通 过 利 | 
FE EH 


EZ LUI 


个 非 封 闭 系统 ， 其 理论 模型 即使 
电离 层 经 验 模 型 虽然 


的 电离 层 结构 。 这 三 个 方向 的 研究 相辅相成 : 
由 于 电离 层 经 验 模 型 精度 较 高 可 作为 同化 模型 的 初始 场 ， 而 电离 层 型 


\ 有 掩 星 电离 层 探测 功能 
不 多 ， 期 待 更 多 的 低 轨 卫星 星座 建成 并 搭载 电离 层 掩 星 天 线 ， 


以 构建 更 加 完善 


t 有 极 高 精度 的 初始 状态 ， 


也 会 随 着 


顾及 了 真实 观测 数据 的 空间 


特征 与 时 间 周 期 


规律 ， 甚 至 可 以 重 现 韦 德 海 现象 ， 但 其 本 质 只 是 一 个 统计 模型 ， 在 极端 空间 


| 环境 依然 有 较 大 


尔 曼 滤波 的 动力 学 部 分 。 以 前 电离 层 观测 数据 主要 来 自 电离 层 垂 测 仪 以 及 


引入 GPS 电离 层 观 测 数据 提高 


然 比 较 少 ， 难 以 实现 磁 暴 或 者 电离 层 ; 


通过 实测 数据 ， 建 立 上 
了 温度 、 离 子 温度 ， 离 子 组 成 等 ， 


电离 层 在 时 间 与 空间 
舌 跃 期 间 的 海洋 区 域 异常 捕获 。 

EE 离 层 气 候 / 经 验 模 型 ， 一 个 基于 数据 的 模型 来 
以 避免 理论 模型 所 产生 的 不 确定 性 (通常 


刁 观 测 数 据 对 电离 层 物理 模式 或 气候 模型 进行 更 新 ， 可 获得 更 加 
以 卡尔 曼 滤波 为 同化 模型 核心 模块 ， 
! 论 模型 可 将 其 作为 卡 


非 相干 散射 雷达 ， 


上 的 探测 能 力 ; 但 电离 层 观测 数据 依 


述 电子 密度 、 电 


有 更 好 的 精 


度 ， 可 以 作为 同化 或 者 层 析 算 法 的 背景 场 )。 目 前 的 电离 层 和 等 离子 体 层 经 验 模型 主要 有 


Nequick" ™, IRI (international reference ionosphere), GCP 


M (global core plasmasphere 


model)", IRI_Plas”” 等 模型 。 电 高 层 F2 层 顶部 区 域 空间 环境 复杂 且 探 测 数据 少 ， 传 统 


模型 存在 较 大 误差 。GNSS 电离 层 掩 星 观测 数据 通过 反 演 可 以 得 到 低 轨 ] 


电子 密度 剖 线 ， 可 ) 


基于 函数 基 或 像素 基 ， 


来 改善 电离 层 顶 部 区 域 经 验 模型 。 
通过 层 析 算法 一 、 卡 尔 曙 滤波” /三维 或 四 维 变 分 /统计 优 


化 ( 解 形式 具有 等 价 性 ””) 等 方式 对 物理 模式 或 气候 模型 进行 同化 更 新 ， 


上 息 ， 但 是 它 并 不 能 提供 电离 层 下 一 刻 时 间 的 误差 评估 ， 也 难以 保 说 


为 了 获得 连续 的 电离 层 三 维 结构 ， 对 电离 层 层 析 结果 进行 平滑 ， 则 会 丧失 
构 。 通 过 引进 可 应 用 于 气象 学 、 
的 误差 评价 ， 以 及 空间 结构 的 连续 性 方面 会 


及 计算 能 


Ep 星 轨 道 高 度 以 下 的 


可 以 获得 精细 准 


确 的 电离 层 结构 。 层 析 算 法 可 以 基于 背景 场 使 用 极 小 的 内 存 与 计算 时 间 获 得 三 维 电离 层 信 


电子 密度 空间 的 连续 性 ; 


已 离 层 局 部 精细 结 
海洋 地 质 学 以 及 地 震 学 的 同化 算法 ， 在 物理 机 制 、 时 间 演 化 
更 好 的 结果 ， 但 同化 模型 对 
RARR. Huang 提出 了 一 种 先进 的 空间 滤波 方案 ， 以 避免 


[于 计算 机 的 内 存 
3D-Var 寻找 背景 


场 误差 协 方差 逆 和 矩阵 ”， 极 大 缩短 了 计算 机 计算 时 间 。Schunk 等 人 使 用 基于 物理 学 的 电离 


ChinaXiv& 'ERBTII 


1 期 NTs, GS: 基于 多 源 GNSS 观测 数据 的 三 维 电 离 层 研究 现状 及 发 展 65 


层 -等 离子 层 模型 和 卡尔 曼 滤 波 器 作为 融合 各 种 近 实 时 测量 资料 的 基础 ， 开 发 了 全 球 电离 层 
同化 模型 (global assimilation of ionosphere measurements, GAIM); GAIM 模型 缺少 GNSS 
多 星座 数据 的 支持 ， 电 离 层 观测 资料 的 时 空 分 布 较 差 ”™”™。Courtier 55 AP" 和 Wang 等 
人 ”将 四 维 变 分 用 于 随时 间 变 化 的 观测 数据 ， 利 用 增 量 近似 方法 为 四 维 变 分 的 成 本 与 收益 
权衡 提供 了 灵活 性 。Yue 等 人 ”将 地 基 及 空 基 GNSS 电离 层 观测 数据 通过 集合 卡尔 曼 滤波 
同化 模型 进行 了 2002 年 到 2011 年 三 维 电离 层 重 分 析 ， 但 是 并 没有 使 用 四 维 变 分 及 电离 层 物 
REUS. Fu eA" 提出 电离 层 同 化 需要 考虑 接收 机 与 发 射 机 的 硬件 延迟 ， 以 及 利用 双 层 分 
步 同化 方式 提高 电离 层 同 化 模型 精度 。 在 多 源 及 多 模 GNSS 电离 层 观 测 数据 的 增加 背景 下 ， 
由 于 多 源 观测 的 多 样 化 数据 特征 及 限于 计算 机 内 存 与 算 力 约束 ， 对 电离 层 同化 模型 构建 过 
程 有 了 新 的 要 求 与 发 展 。 

第 回 章 介绍 了 多 模 地 基 及 星 载 GNSS 电离 层 观 测 数 据 发 展现 状 。 第 团 章 讨论 了 空 
GNSS 电离 层 观 测 数据 在 电离 层 顶部 区 域 经 验 模型 改进 方面 的 进展 。 oe 
模 GNSS 电离 层 观测 数据 的 增加 背景 下 电离 层 同 化 模型 构建 过 程 中 新 的 要 求 与 发 展 。 第 回 
章 提出 自 适 应 格 网 与 机 器 学 习 可 成 为 电离 层 同 化 算法 方面 的 新 研究 方向 ， 精 细 三 维 电离 层 
结构 通过 利用 理论 模型 /经 验 模型 与 多 源 电 离 层 观 测 数据 进行 融合 有 了 极 大 提高 ， 并 被 付 诸 
诸多 方面 应 用 与 研究 ， 以 促进 三 维 电离 层 研究 。 


2 多 源 GNSS 电离 层 观 测 数据 


电离 层 数据 依照 数据 特征 或 探测 方式 可 以 分 为 以 下 三 种 。 一 种 是 太阳 活动 、 磁 暴 情 
况 、 地 磁 变 化 等 多 种 指 征 性 数据 ， 如 地 磁 指 数 ， 地 磁场 扰动 幅度 /地 磁 行星 指数 ， 太 阳 黑 子 
数 等 一 一 它们 只 是 对 电离 层 及 地 球 空间 环境 做 出 宏观 的 描述 。 另 一 种 是 通过 直接 探测 手段 
获得 的 电离 层 数 据 ， 如 来 自 于 电离 层 垂直 探测 仪 (ionosonde)、 非 相干 散射 雷达 (incoherent 
scatter radar, ISR) 等 手段 的 数据 。 这些 方式 覆盖 范围 小 、 探 测 成 本 高 、 探 测 时 间 短 且 连 
续 性 低 ， 但 获得 的 数据 精度 较 高 ， 可 以 作为 电离 层 模 型 构建 及 验证 的 信息 。 最 后 一 种 是 
tH GNSS 技术 发 展 以 来 ， 通 过 从 微波 穿 过 电离 层 得 到 的 GNSS 载波 相位 及 伪 距 中 解析 出 
电离 层 影响 ， 进 而 对 电离 层 电子 密度 分 布 进行 建 模 获 得 的 相关 电离 层 产品 ， 如 国际 GNSS 
服务 (International GNSS Service, IGS) 的 两 维 电 离 层 图 GIM (global ionosphere model) 产 
m- COSMIC (constellation observing system for meteorology, ionosphere, and climate) 系 
统 通 过 CDAAC/UCAR (the COSMIC Data Analysis and Archive Center (CDAAC) of the 
University Corporation for Atmospheric Research) 组 织 提 供 的 电离 层 电 子 密度 廓 线 ion Prf 
数据 。 以 下 将 介绍 基于 GNSS 技术 间接 探测 获得 的 电离 层 数据 特征 、 
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基于 地 基 GNSS 台 站 、 低 轨 卫 星 和 GNSS 导航 卫星 ， 图 四 展示 了 空间 科学 在 大 气 / 海 
洋 /电离 层 等 方面 的 观测 系统 。 以 电离 层 为 例 ， 地 基 GNSS 卫星 系统 及 地 面 监测 台 站 构成 
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了 地 基 GNSS 电离 层 观 测 系统 ， 低 轨 卫 星 /GNSS 卫星 构成 了 高 层 电 离 层 观 测 及 掩 星 电 离 层 
观测 系统 ， 低 轨 卫 星 间 通过 相互 发 射 与 接收 信号 形成 低 轨 卫星 间 的 电离 层 监 测 系 统 ， 例 如 
COSMIC/COSMIC-2 (分 别 发 射 于 2006 和 2019 年 );， 以 Jason-1/Jason-2/Jason-3 (分 别 发 
射 于 2001, 2008, 2016 年 ) 为 例 ， 通 过 海洋 垂直 观测 形成 了 海洋 电离 层 观 测 系统 。 


—— Ground-based GNSS 
== GNSS-LEOOccult 
一 一 GNSS-LEO Overhead 
一 一 LEO-LEO CrossLink 
GNSS/LEO Reflection 


Jason-Sea Observation 


1 ”现在 及 将 来 中 性 大 气 /海洋 /电离 层 数 据 观 测 系统 


2.1 地基 GNSS 电离 层 观 测 

AO Ata, AINA) GNSS 卫星 与 地 面 台 站 间 信 和 号 传递 过 程 中 形成 标 色 线 指示 的 路 
径 ， 它 跨越 了 平流 层 、 对 流 层 、 电 离 层 以 及 等 离子 体 层 。 卫 星 信 号 在 不 同 圈 层 穿 过 时 ， 由 
于 不 同 圈 层 的 物理 特性 对 微波 信号 产生 各 种 附加 影响 ， 虽 然 给 GNSS 定位 带 来 了 许多 误差 ， 
但 也 为 GNSS 数据 反 演 各 个 圈 层 的 物理 特性 提供 了 条 件 。 
现在 已 有 完善 的 GNSS 卫星 星座 ， 包 括 GPS, GLONASS (Global’naya sputnikovaya 
navigatsionnaya sistema), Galileo, Beidou 等 导航 系统 。 各 系统 卫星 具有 不 同 轨道 设计 方式 
及 特征 ， 使 得 混合 系统 可 以 保证 地 面 上 每 个 台 站 至 少 有 十 几 个 卫星 可 见 。 以 地 基 GNSS / 
IGS 观测 系统 为 例 ，30 s 采样 率 的 观测 数据 每 小 时 可 以 得 到 数 十 万 的 观测 数据 。 图 回 展 
示 了 2019 年 IGS 以 及 中 国 陆 态 网 ， 澳 大 利 亚 GeoScience、 新 西 兰 GeoGoing、 美 国 /日 本 
GeoNet 的 区 域 监测 网 的 各 站 位 置 分 布 。 

利用 全 球 地 基 GNSS 观测 ， 人 们 可 以 获得 海量 的 全 天 候 、 全 球 分 布 的 电离 层 观 测 数据 。 
如 图 加 地基 GNSS 观测 台 站 主要 分 布 在 大 陆 上 ， 少 量 建 在 海岛 上 ， 叶 致 在 海洋 区 域 的 地 基 
GNSS 电离 层 观测 数据 十 分 稀少 。 如 果 要 建立 全 球 电离 层 模型 ， 需 要 补充 海洋 区 域 的 观测 
数据 。 
2.2” 空 基 GNSS 电离 层 观 测 

空 基 GNSS 电离 层 观 测 数据 主要 分 为 三 种 : (1) 用 于 卫星 定 轨 过 程 中 的 定位 数据 可 以 获 
得 卫星 到 导航 星之 间 斜 向 电子 总 量 ; (2) 卫星 到 导航 星之 间 发 生 掩 星 时 可 以 得 到 电离 层 掩 星 
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2 ”2019 年 全 球 IGS 站 及 区 域 监测 网 分 布 


观测 数据 ; (3) 测 高 类 卫星 可 以 获得 卫星 高 度 到 地 球 表面 之 间 的 垂直 电子 总 量 。 从 电离 层 观 
测 的 角度 来 看 ， 基 于 电离 层 临 边 探测 的 捧 星 数据 具有 高 精度 、 高 垂直 分 辨 率 、 完 整 的 全 球 
覆盖 以 及 无 电离 层 偏差 等 优势 。2006 年 中 国 台 湾 和 美国 联合 发 射 的 6 颗 低 轨 卫 星星 座 组 成 
的 COSMIC 系统 为 电离 层 研究 提供 了 大 量 的 掩 星 观测 数据 。 当 然 近年 来 由 于 COSMIC Æ 
座 运 行 寿命 已 尽 ， 只 能 提供 很 少 的 观测 数据 ， 但 是 采用 GPS 和 GLONASS 双 模 观测 的 6 颗 
低 轨 卫星 星座 COSMIC-2 已 于 2019 年 发 射 ， 将 提供 更 多 的 空 基 电离 层 观 测 数据 。 如 图 加 所 
示 ，CDAAC/UCAR 近 实 时 处 理 和 存档 10 多 个 掩 星 任务 数据 ， 提 供 了 进行 全 球 三 维 电 离 层 
电子 密度 重 构 及 重 分 析 的 一 个 机 遇 。 图 田 使 用 表 四 中 的 参数 模拟 了 2008 年 01 月 10 日 全 天 
328 个 IGS 观测 站 对 GPS 和 GLONASS 两 个 星座 系统 的 电离 层 地 基 观 测 数据 的 电离 层 穿 刺 
点 (ionospheric piercing points, IPP), LAR COSMIC 六 个 卫星 对 GPS 星座 的 电离 层 掩 星 观 
测 数据 的 电离 层 切 点 (the target point of the occultation signal path, OccTarget) 的 分 布 情 
况 。 从 图 田中 可 知 ， 电 离 层 掩 星 观测 数据 可 以 极 大 地 改善 电离 层 地 基 观 测 数据 在 海洋 区 域 观 
测 匮 乏 问 题 ， 然而 相对 于 众多 电离 层 地 基 观 测 数据 ， 仅 由 COSMIC 低 轨 卫星 星座 产生 的 电 
离 层 掩 星 观测 数据 较 少 ， 导 致 海陆 观测 资料 分 布依 然 不 均衡 。 


= 


表 1 地 基 / 星 载 电离 层 观 测 
模拟 系统 地 基 星 载 
导航 星座 及 卫星 数 ”GPS +GLO(56) GPS (32) 
接收 机 载体 及 数目 IGS (256) COSMIC (6) 


截至 高 度 角 /(? ) 15 = 
采样 率 /s 300 5 
观测 数目 /h 约 10000 约 10000 


近 5 年 内 ，OneWeb、SpaceX 和 Boeing、Samsung、 航 天 科 工 集团 及 航天 科技 集团 等 ， 
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Jason-2 


ccs FIT ~~ 


T T TT] T 1 E 1 LIT T DT | TITY TT 


UCAR/CDAAC 日 累计 掩 星 事 件数 目 


探测 的 数 


2001 2002 2003 2004 2005 2006 2007 2008 2009 2010 2011 2012 


年 份 


图 3 2001—2012 年 间 携带 掩 星 有 效 载荷 的 卫星 任务 掩 星 数据 


* IPP 
* Occ Target 


图 4 ”全 天 电离 层 观测 数据 分 布 情况 


202306.00339v1 


chinaXiv 


ChinaXiv 合 作 期 


1 期 付 乃 锋 ， 等 ， 基 于 多 源 GNSS 观测 数据 的 三 维 电 离 层 研究 现状 及 发 展 69 


为 了 基于 太空 向 全 球 提供 无 缝 旦 稳定 的 互联 网 服务 ， 各 自 先 后 设计 、 部 署 、 并 计划 发 射 数 十 
至 上 万 颗 不 等 的 商业 低 轨 卫星 星座 。Zhang 和 Ma ”整理 了 各 卫星 星座 设计 参数 ， 如 表 回 所 
示 。 而 Reid 等 人 ”从 整个 系统 的 架构 ， 包 括 星座 几何 图 形 条 件 、 空 间 信号 测 距 误 差 、 星 载 
原子 钟 性 能 和 定 轨 方 法 等 方面 ， 全 面 探索 其 扩展 成 为 导航 增强 星座 的 可 能 性 ， 结 果 令 人 振 
奋 。 电 离 层 撼 星 载荷 未 来 小 型 化 与 低 成 本 化 ， 将 其 搭载 在 不 同 导 航 增强 星座 上 ， 可 以 形成 分 
布 更 均匀 ， 数 目 更 多 的 电离 层 掩 星 事 件 ， 进 而 获得 更 多 电离 层 掩 星 观测 数据 用 于 电离 层 三 维 
重 构 。 


xo ”部 分 已 部 署 或 提出 的 商用 低 轨 星 座 汇 总 


星座 卫星 数 “高度 ”倾角 ”建成 时 间 国家 主要 业务 
/km  /(^) /年 份 
Iridium 66 780 86.4 1998 美国 语音 十 STL 
Globalstar 48 1400 52 2000 美国 语音 
Iridium NEXT 75 780 | 86.4 2019 美国 宽带 十 STL 
648 1200 pem 
OneWeb Vues = 88 2027 美国 宽带 


1600 1150 53 
1600 1110 53.8 


SpaceX Starlink 2024 美国 宽带 


7518 340 

1 190 45 
Boeing 612 1200 55 - 美国 宽带 

1 155 88 
LeoSat 108 1400 — 2020 美国 宽带 

72 1000 99.5 要 M 

Telesat 45 lodge d 2022 加 拿 大 宽带 
Kepler 
Comunications 140 -一 一 2022 加 拿 大 物 联 网 
Astrocast 64 600 一 2021 瑞士 物 联 网 
Yaliny 135 600 一 一 dU 宽带 
Astrome 150 1400 — 2020 印度 宽带 
Samsung 4600 1400 — — 韩国 宽带 
Hongyan 1100 E SUR rp ER 宽带 十 导航 增强 
Hongyun 156 1000 — 2022 中 国 ”宽带 二 导航 增强 
CentiSpace 120 700 =m 2021 中 国 ”宽带 上 + 导航 增强 
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3 多 源 观 测 GNSS 数据 应 用 于 电离 层 经 验 模 型 


从 SAMI2/3 模型 中 四 到 NCAR/TIEGCM 模型 下 ， 再 到 覆盖 从 地 面 到 磁 层 的 整 
体 模型 (ground-to-topside model of atmosphere and ionosphere for aeronomy, GAIA)"", 
及 已 /正在 扩展 的 包括 带 有 热 层 和 电离 层 扩 展 的 整个 大 气 的 气候 模型 (the thermosphere and 


人 


ionosphere extension of the whole atmosphere community climate model, WACCM-X 
电离 层 理论 模型 获得 了 长 足 的 发 展 。 以 GATA 模型 为 例 ，GAIA 是 三 个 独立 开发 的 模型 的 组 
合 ， 包 括 整个 大 气 层 (从 对 流 层 到 热 层 )、 一 般 环 流 模 型 (GCM)、 电 离 层 模 型 和 电动 力学 模 
Al. GAIA 的 核心 是 一 个 “耦合 器 ”模块 ， 用 于 管理 三 种 模型 之 间 的 差异 。 人 
大 气 层 的 模型 (如 TIE-SGCM), RE GAIA 模拟 了 从 地 球 表 面 到 热 圈 的 中 性 大 气 ”。 相 对 而 
言 ， 电 离 层 理论 模型 作为 同化 模型 的 初始 场 并 没有 电离 层 经 验 模 型 精度 高 ， 往 往 将 其 作为 卡 
尔 曼 滤波 的 动力 学 部 分 ， 所 以 利用 多 源 数据 改进 电离 层 经 验 模型 是 当前 的 研究 热点 之 一 。 由 
于 电离 层 F2 层 顶部 区 域 空间 环境 复杂 且 探 测 数据 少 ，IRI 模型 、GCPM 模型 以 及 IRI Plas 
模型 在 电离 层 顶部 建 模 做 了 大 量 工作 。 这 里 主要 介绍 IRI 模 型 及 基于 掩 星 资 料 改进 TRI 模型 
的 研究 进展 。 

国际 参考 电离 层 TRI 基于 对 地 面 和 空间 电离 层 等 离子 体 的 大 多 数 可 用 和 可 靠 的 观测 结 
果 ， 描 述 了 高 度 60—2000 km 范围 内 的 电子 密度 、 电 子 温 度 、 离 子 温度 、 离 子 组 成 以 及 其 
他 参数 的 月 平均 值 。 相 比 IRI2012， 最 新 版 本 的 IRI2016 有 两 个 新 的 关于 F2 层 峰 值 高 度 选 
项 ， 以 后 不 再 依赖 与 传播 因 了 M(3000) 的 关系 来 描述 F2 层 峰值 高 度 ; 个 是 基于 电离 层 测 
高 仪 数 据 ”， 被 推荐 以 获得 更 加 准确 的 电离 层 的 电子 密度 分 布 场 ”， 另 一 个 基于 掩 星 观 测 数 
m, 

相 比 于 地 基 GNSS HA, HS Ji d EE ULLA HES ORE DE 
来 ， 许 多 研究 人 员 做 了 大 量 研究 以 结合 COSMIC HH BS Het E LU dE BS TP 
模型 。 采用 经 验 正 交 函 数 (empirical orthogonal function, EOF) 分 析 方 法 ， 处 理 2007 一 2011 
年 间 COSMIC 电离 层 掩 星 数据 并 重建 了 电离 层 标高 模型 ， 其 中 地 磁 纬 度 (87.5°S—87.5°N) 
的 空间 分 ee SHEN 2 h。 通 过 EOF 分 析 ， 发 现 标高 的 特征 主要 体现 在 地 磁 
纬度 、 年 度 、 季 节 和 昼夜 变化 中 ， 如 图 加 所 示 ， 通 过 将 比例 因子 q 作为 附加 约束 改进 国际 
55 HR E 依 I9001 问 的 顶部 模型 ， 这 一 部 分 曲线 拟 合 精度 得 到 极 大 提高 ”， 并 将 这 一 模 
型 辅助 Abel 反 演 以 验证 其 精度 ”。 另 一 方面 ， 在 2018 年 通过 COSMIC 的 顶部 总 电子 密度 
观测 及 映射 函数 得 到 垂直 电子 密度 总 量 ， 再 利用 垂直 电子 密度 总 量 与 来 自 COSMIC 电离 层 
星 数 据 中 的 卫星 轨道 高 度 处 的 电子 密度 和 顶部 区 域 标高 关系 ， 建 立 了 800 km 以 上 区 域 的 
电离 层 电子 密度 指数 衰减 模型 ， 以 用 于 电离 层 顶部 及 等 离子 体 层 影响 的 评价 ， 或 误差 的 消 
ps, 


基于 COSMIC 电离 层 掩 星 电子 密度 曲线 数据 建立 的 电离 层 标高 模型 作为 附加 约束 改善 
SIRT 的 顶部 模型 ， 但 恒定 的 标高 不 适合 描述 电离 层 电子 密度 的 物理 特性 ， 尤 其 在 等 离子 体 
层 会 低估 电子 密 ". 因此 ， 利 用 Vary-Chap 变 标 高 (Vary-Chap scale height, VCSH) 技术 有 
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900 900 
— — COSMIC —— COSMIC — COSMIC 
RAUM —— — — non cor 
nm 800 cor 800 = 
700 700 
600 600 
g 500 500 
m 
E 400 400 
300FHSC 126 km 300 
200 200} 
100 100 
> Tt 6 $ oxi ? 1 2 axip 0— 0 —8 10 10! 
包子 密度 /cms ET EE em? SE RE / em? 
a) b) c) 


iE: a) HX 2007.030, HAX 10.1958, EX 24.97, ZEX 105.40; b) 日 期 为 2007.030， 世 界 时 
为 10.2158， 纬 度 为 34.75， 经 度 为 95.66; c) 日 期 为 2007.030， 世 界 时 为 10.2800， 纬 度 为 40.07， 经 度 为 
129.85. “HSC” KI F2 层 顶部 标高 ,“non cor” 表 示 q 因子 约束 校正 前 ,“cor” 表 示 q 因子 约束 校正 后 。 


图 5 ”比例 高 度 约束 引入 IRI 前 后 ， 三 个 反 演 结果 与 检索 对 应 的 COSMIC 模拟 轮廓 的 比较 结果 


望 进一步 改善 电离 层 顶部 电子 密度 模型 。Vary-Chap 变 标高 H(h) 于 1969 年 首先 被 引入 到 


一 般 a-Champman 电子 密度 分 布 : 


Ne(h) = Ny (de) exp [5 - sc epto3)] ， 0) 
^ dh 
之 二 n. H(h) , (2) 


Hp, 2A BEES Nm, hm TIRE F2 层 峰 值 密度 及 高 度 ，N(h) RATER h SO 
HP; H(h), Hm 是 高 度 hh 及 F2 层 峰 值 高 度 处 的 标高 ， 即 Hm = H (hm) H(h) BRA 
ARE Hm 的 函数 。 利 用 2008 *E 1H 1 H £8 2013 4 12 A 31 Hf COSMIC Hà Tx E Ze 
数据 ， 得 到 两 个 由 指数 形式 的 权重 控制 的 线性 az 形状 与 抛物 线 ba? 形状 模式 组 合 的 形状 函 
数 ， 拟 合 代 表 较 低 和 较 高 高 度 的 Vary-Chap 标高 廊 线 : 


y = 0 一 一 一 一 F pre, (3) 
Joh, y= In(H(h)) —In(Hy), z = “7 是 中 间 变量 ， 用 以 建立 标高 与 高 度 之 间 的 关系 ; 
a 与 5 分 别 是 线性 与 抛物 线 模式 形状 的 拨 合 参数 。 当 使 用 具有 权重 的 线性 模式 ar 和 抛物 线 
模式 ba? 来 表示 VCSH 轮廓 的 形状 时 ， 拟 合 的 VCSH 和 重建 的 电子 密度 轮 廊 与 原始 观测 值 
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非常 吻合 ; 在 所 有 研究 的 年 份 中 ， 大 多 数 拟 合 残 差 都 小 于 10%。 参 数 a 的 范围 为 1.0 一 2.0 
(每 月 平均 值 )， 与 地 磁场 线 高 度 相 关 ， 并 且 在 白天 和 高 活动 年 位 于 低 纬度 和 中 纬度 区 域 较 
K; 而 参数 5 与 转换 高 度 (两 种 形状 模式 的 交点 ) 的 变化 规律 与 磁 亦 道 异常 (the equatorial 
ionization anomaly, EIA) 极为 相似 ”。a 和 尹 参 数 化 拟 合 函数 可 以 将 VCSH 轮廓 以 合理 的 
精度 简单 表示 ， 可 以 结合 F2 层 的 峰值 电子 密度 、 峰 值 高 度 和 标 度 高 度 在 任何 给 定 的 时 间 和 
位 置 重建 电子 密度 分 布 ， 经 过 仔细 验证 ， 它 可 以 作为 另 一 个 topside 选项 引入 IRI 模型 ”。 


4 ZI GNSS 观测 数据 应 用 于 电离 层 同 化 模型 


近年 来 ， 基 于 代数 重建 技术 ， 电 离 层 层 析 成 像 算法 已 广泛 用 于 电离 层 3D 结构 的 重 
建 ”™。 然 而 ， 电 离 层 的 观测 数据 空间 分 布 差 ， 层 析 成 像 结 果 缺 乏 空 间 连 续 性 。Lorenc 
mbox” 采 用 卡尔 曼 滤波 方法 提出 了 三 维 变 分 同化 作为 现代 数据 同化 的 经 典 方法 。 图 日 展示 
了 同化 后 的 电离 层 三 维 模型 在 经 度 、 纬 度 及 高 度 上 的 切片 ， 它 更 细致 地 描述 出 电离 层 三 维 结 
构 特 性 ， 特 别 是 电离 层 赤 道 附近 的 不 规则 形状 。 


Elo 电离 层 电 子 密度 场 沿 经 度 、 纬 度 及 高 度 上 的 切片 


与 层 析 成 像 算法 相似 ， 电 离 层 同化 算法 利用 观测 值 来 更 新 背景 场 : 基于 电离 层 电子 密 
度 背 景 场 及 观测 数据 各 自 的 协 方差 ， 构造 和 求解 与 观测 值 和 背景 场 有 关 的 价值 函数 极 小 值 ， 
然后 得 到 同化 结果 。 该 方法 可 以 保持 电离 层 的 光滑 物理 结构 ， 并 且 更 接近 真实 电离 层 分 布 包 
电子 密度 场 。Huang 提出 的 空间 滤波 方案 避免 了 3D-Var 寻找 背景 场 误差 协 方差 着 矩阵 ” ， 
但 是 多 源 GNSS 电离 层 观测 数据 的 增多 ， 依 然 对 计算 机 的 内 存 与 计算 能 力 有 极 高 要 求 。 采 
用 逐步 双 层 同化 方法 来 减少 电离 层 项 部 区 域 对 电离 层 底部 区 域 同化 的 影响 ”: 同化 低 轨 卫 
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星 定位 观测 获得 斜 向 TEC 信息 ， 得 到 电离 层 800 km 以 上 区 域 电子 密度 分 布 模型 ， 再 同化 
掩 星 观测 数据 与 扣除 800 km 以 上 电离 层 的 地 面 GNSS 电离 层 观测 ， 获 得 800 km 以 下 的 电 
离 层 电子 密度 分 布 模型 ， 以 得 到 更 准确 的 F 层 电离 层 结构 。 这 种 方式 可 以 减少 分 步 同化 中 
未 知 量 数目 ， 进 而 节省 计算 时 间 与 计算 机 内 存 。 
可 以 看 到 ， 多 源 及 多 模 GNSS 电离 层 观测 数据 的 引入 ， 使 得 电离 层 同化 模型 有 了 新 的 
要 求 与 发 展 : (1) 为 保证 多 源 及 多 模 GNSS 电离 层 观测 数据 特征 及 精度 ， 同 化 模型 中 各 种 误 
差 源 需 要 更 细致 的 分 析 ， 并 寻求 方法 加 以 消除 ，(2) 需要 对 多 源 及 多 模 GNSS 电离 层 观测 数 
据 在 同化 模型 中 的 改进 效果 进行 评价 ， 由 于 观测 数据 增多 要 求 算法 设计 需 适 应 计算 效率 。 
4.1 ”同化 算法 中 的 误差 分 析 
在 电离 层 同化 过 程 中 ， 背 景 场 误差 分 布 的 假设 、 观 测 数据 的 穿刺 路 径 分 布 、 观 测 数据 的 
权重 ， 以 及 对 和 抢 阵 解 进行 稀 疏 优化 等 都 可 能 导致 同化 结果 出 现 误 差 。 对 误差 源 进行 分 类 并 有 
效 地 改善 这 些 误 差 是 提高 同化 结果 精度 的 重要 手段 ， 可 将 电离 层 误 差 源 分 为 算法 误差 与 模 
型 假设 误差 。 
= 虽然 算法 误差 难以 扣除 ， 通 过 模拟 已 验证 ， 经 验 性 调整 电离 层 背景 场 误差 及 相关 系 
Ey 数 、 背 景 场 权重 与 各 不 同 观测 误差 ， 同 化 结果 误差 可 降低 到 真实 场 的 20%， 而 整体 偏差 小 
于 5%。 然 而 几 个 来 源 于 同化 过 程 假设 的 误差 ， 使 得 模拟 观测 数据 与 固定 背景 场 误差 具有 大 
的 系统 偏差 ， 但 已 验证 可 以 通过 算法 进行 改进 。 如 图 四 所 示 ， 电 离 层 顶层 及 等 离子 体 层 影 
响 、 电 离 层 时 间 变 化 影响 、 格 网 代表 性 误差 等 可 以 通过 有 效 的 方法 减弱 ， 都 属于 模型 假设 误 
差 。 以 电离 层 顶层 及 等 离子 体 层 影响 为 例 ， 这 一 部 分 对 GPS 的 TEC 的 贡献 随 一 天 中 的 时 
间 变 化 ， 从 昼夜 平分 的 最 低 (25 12%) 到 冬天 夜间 的 最 大 ( 约 60%)”。 而 以 GCPM 模型 进 
行 模拟 统计 后 ， 发 现 所 有 模型 假设 误差 一 般 情况 下 可 能 达到 10%~20%， 与 诸多 学 者 结论 
致 ™ 修 ， 扣 除 电 离 层 顶 层 及 等 离子 层 影响 后 ， 则 减弱 到 8% 左右 。 


* Satellite02 * Satellite03 


* Satellite01 


ik: “STECXX” # an # lal rm dm ROUGH (XX 表示 编号 ),“SatelliteXX” 表 示 GNSS 导航 星 ， 
“ReceiverXX ”表示 地 面 观测 站 。 


7 ”电离 层 顶层 及 等 离子 层 影 响 、 电 离 层 时 间 变 化 影响 、 格 网 代表 性 误差 对 同化 过 程 的 综合 作用 示意 图 
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对 于 电离 层 顶部 及 等 离子 体 层 影 响 ， 掩 星 数 据 可 以 用 于 顶部 区 域 来 获得 更 准确 的 顶层 
电子 密度 衰减 函数 的 基数 ”， 而 COSMIC 项 目 中 的 LEO 卫星 轨道 高 度 约 为 800 km， 所 以 
接近 卫星 轨道 高 度 的 校正 电子 总 量 可 以 作为 电离 层 同 化 的 约束 ， 从 而 实现 顶部 和 底部 区 域 
同化 模型 之 间 的 一 致 性 和 连续 性 。 在 处 理 观 测 数据 时 ， 选 择 顶 部 区 域 和 底部 区 域 用 于 逐步 
同化 。 首 先 ， 将 顶部 区 域 的 经 验 模型 (u IRL Plas”) 用 作 背 景 场 ， 并 通过 提供 podTec 数 
据 获得 顶部 区 域 的 同化 模型 。 然 后 ， 从 地 面 观测 数据 中 减 去 顶部 区 域 贡献 ， 其 与 800 km 以 
下 的 电离 层 经 验 模型 作为 底部 区 域 的 背景 场 进行 同化 后 ， 获 得 底部 区 域 的 高 精度 电离 层 三 
维 电子 密度 场 ”。 其 同化 结果 表明 ， 在 减 去 电离 层 顶 部 的 影响 后 ， 平 静 期 的 电离 层 偏差 从 
1.645TECu 改善 到 1.464 TECu; 当 电 离 层 处 于 活跃 期 时 ， 标 准 偏差 从 4.408TECu 降低 到 
3.536 TECu. 

虽然 电离 层 在 同化 时 间 段 内 变化 明显 ， 但 一 般 在 电离 层 反 演 中 假定 在 1h 或 2h 内 是 固 
定 不 变 的 ， 这 导致 电离 层 反 演 过 程 中 因为 时 间 变 化 引起 的 误差 增 大 了 电离 层 同 化 结果 的 误 
差 。 减 少时 间 窗 口 可 以 提高 同化 结果 的 准确 性 ， 还 可 以 减少 每 次 同化 的 观测 次 数 ， 从 而 减少 
每 种 同化 算法 的 执行 时 间 ”。 但 是 ， 减 少 观 测 值 可 能 会 影响 其 在 电离 层 网 格 中 分 布 的 均匀 
性 。 因 此 ， 需 要 找到 一 种 减轻 同化 窗口 时 间 内 电离 层 变化 影响 的 方法 。 忽 略 电离 层 赤道 异常 
情况 下 在 地 磁 坐 标 下 随 太 阳 直 射 点 移动 ， 电 离 层 分 布 具 有 一 定 的 时 间 平 移 性 ， 所 以 通过 射线 
时 间 与 电离 层 背 景 场 时 间 的 较 差 ， 通 过 公式 进行 地 磁 经 度 平 移 可 获得 校正 后 射线 坐标 ; 


Pmag = Pmag T (M J Dbackground T MJ Dray) x 360° , (4) 


HF, Ónag mag 分 别 是 校正 前 后 的 地 磁 经 度 ，M .JD@sy，M JDbaccground 分 别 是 观测 数据 
和 电离 层 背 景 场 时 间 对 应 的 儒 略 日 期 。 经 过 统计 相对 精度 发 现 ， 通 过 地 磁 坐 标 校正 时 间 偏 差 
的 方法 校正 不 同时 刻 的 倾斜 TEC， 其 结果 在 一 天 之 内 每 小 时 都 产生 稳定 的 积极 影响 ， 电 离 
层 时 间 变 化 的 影响 为 9.27%， 经 过 地 磁 坐 标 时 间 校 正 后 ， 降 低 到 5.62% (4) 0.3 TECu). 
在 电离 层 格 网 化 中 ， 经 典 路 径 妃 踪 方 式 通过 射线 穿 过 每 个 格 网 的 距离 形成 了 整个 射线 
路 径 ， 然 后 通过 每 个 格 网 截 距 及 其 对 应 格 网 中 心 的 电子 密度 来 表示 每 段 路 径 贡 献 的 电子 量 。 
然而 当 穿 刺 路 径 处 于 格 网 一 隅 ， 或 是 四 角 ， 或 是 极为 接近 另 一 个 格 网 时 ， 再 次 以 当前 格 网 中 
心 的 电子 密度 来 表示 此 段 路 径 上 的 平均 密度 并 不 合适 ， 而 穿刺 截 距 的 中 点 p 所 在 位 置 的 电子 
密度 可 以 更 精确 地 代表 此 段 路 径 上 的 密度 。 点 2 的 密度 可 以 通过 其 周边 格 网 中 心 密度 加 权 表 
示 ， 所 以 我 们 可 以 得 到 以 下 关系 式 : 


N Fsegment NE, x ` wiN Erz, , (5) 


FEA, NE scgment; N Ep, N Ez, 分 别 是 截 距 平均 密度 、 截 距 中 点 密度 及 周围 网 格 的 密度 ，wwi 
是 每 个 密度 的 权重 。 经 过 统计 双 线 性 系数 加 权 表 示 格 网 截 距 的 密度 后 ， 发 现 原来 电离 层 网 格 
表示 的 误差 约 为 2.5%， 但 使 用 双 线 性 插值 方法 后 降低 到 0.996. 

图 四 模拟 了 GPS 地 基 观 测 数据 空间 分 布 、 电 离 层 顶部 和 等 离子 体 层 影响 、 电 离 层 时 间 
变化 、 格 网 代表 性 误差 三 种 模型 假设 误差 影响 校正 前 后 同化 GPS 地 基 观 测 数据 获得 的 电离 
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层 密 度 分 布 模型 精度 。 由 图 8a) 可 见 ， 地 基 GPS 电离 层 观测 数据 主要 分 布 在 大 陆 及 一 些 岛 
屿 附近 ; 由 图 8b),c) 可 知 ， 忽 略 模型 假设 误差 与 同化 地 基 电 离 层 观测 数据 获得 的 电离 层 密 
度 分 布 模型 相 比 ， 同 化 模型 初始 场 在 全 球 范 围 内 精度 有 明显 提高 ， 在 赤道 附近 及 海洋 区 域 有 
较 大 误差 图 8c),d) 说 明 地 基 GPS 电离 层 观 测 数据 经 过 模型 假设 误差 校正 之 后 ， 同 化 地 基 
电离 层 观 测 数据 获得 的 电离 层 密度 分 布 模型 在 南美 洲 及 赤道 附近 的 误差 大 幅度 降低 了 。 


100 « 


m 


S 
& 
截 距 长 度 平 均 密 度 /km'c 


垂直 TEC 平 均 误 差 /TECu 


B c) 


ii: a) GPS 地 基 观 测 数据 水 平 穿刺 信息 ，b) 初始 场 垂直 TEC 的 平均 误差 ，c) 假设 误差 校正 前 垂直 TEC 的 
平均 误差 ，d) 假设 误差 校正 后 垂直 TEC 的 平均 误差 。 


8 ”日 均 水 平 穿刺 信息 及 垂直 TEC 的 平均 误差 分 布 


4.2 ”多 源 及 多 模 数据 的 使 用 
基于 全 球 导 航 卫 星系 统 对 地 面 站 电离 层 总 电子 含量 的 观测 ，CHAMP (challenging 
minisatellite payload), GRACE (gravity recovery and climate experiment), COSMIC, SAC- 


C (satellite de aplicationes cientificas-B), Metop-A (meteorological operational satellite - A) 
和 TerraSAR-X 卫星 对 电离 层 无 线 电 掩 星 的 观测 ，Jason-1/Jason-2 垂直 电子 总 量 测 量 等 多 
源 数 据 ，2012 年 通过 国际 参考 电离 层 TRI2007 模型 和 卡尔 曼 滤 波 技术 构建 了 2002—2011 年 
全 球 三 维 电 离 层 数 据 同化 模型 ”。 同 时 COSMIC 的 顶部 斜 向 电子 总 量 也 被 用 来 建立 顶部 
垂直 电子 总 量 模型 ， 或 者 通过 同化 和 层 析 的 方式 来 建立 电离 层 顶部 及 等 离子 体 层 电子 密度 
gUu7 7. COSMIC 运行 寿命 将 尽 ， 但 COSMIC-2 星座 计划 (在 500 km 高 度 上 拥有 6 颗 
24° 倾角 卫星 ， 在 800 km 高 度 上 拥有 6 W72 倾角 卫星 ) 可 以 更 好 地 替代 。 通 过 设计 和 模拟 
掩 星 观测 数据 ， 并 基于 经 验 国际 参考 电离 层 模型 和 卡尔 曼 滤 波 器 ， 建 立 了 一 个 高 度 范 围 在 
80—3000 km 的 全 球 电离 层 数据 同化 模型 。 同 化 结果 验证 模拟 的 COSMIC-2 掩 星 观 测 数 据 
补充 了 现 有 的 地 面 GNSS 观测 网 ， 同 时 24° 和 72° 倾角 卫星 可 以 互补 ， 优 化 全 球 电 离 层 建 
WC. RA GPS GLONASS 双 模 观测 的 COSMIC-2 计划 已 于 2019 年 发 射 ， 但 由 于 经 费 
问题 ， 仅 由 6 颗 履 盖 中 低 纬 的 24° 倾角 低 轨 卫星 星座 组 成 ， 仅 为 中 国 台湾 提供 空间 服务 。 基 
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于 表 四 中 的 参数 进行 GNSS 电离 层 观测 数据 模拟 ， 图 回 展 示 了 单一 GPS 星座 地 基 观 测 数据 
与 多 模 地 基 和 掩 星 观 测 数 据 的 空间 分 布 对 比 ， 以 及 数据 应 用 于 同化 的 结果 精度 改进 情况 。 由 
9a),b) 可 知 ， 多 模 地 基 及 掩 星 观 测 数据 在 空间 分 布 要 比 单一 GPS 地 基 电 高层 观测 数据 量 
更 多 ， 禾 盖 范围 更 全 面 ， 由 图 9c),d) 可 知 ， 同 化 单一 GPS 星座 地 基 电 离 层 观测 数据 获得 的 
电离 层 模型 相 比 同化 模型 初始 场 在 全 球 范围 内 精度 有 明显 提高 ， 海 洋 区 域 由 于 观测 相对 较 
少 依然 有 较 大 误差 ， 由 图 9d),e) 可 知 ， 多 模 地 基 及 掩 星 观测 数据 同化 结果 与 GPS 地 基 观 测 
数据 的 同化 结果 相 比 ， 前 者 在 海洋 区 域 的 误差 更 小 ， 但 是 海洋 区 域 依 然 有 较 大 误差 ， 需 要 更 
多 电离 层 GNSS 掩 星 观测 数据 。 


100 s 


m 


心 
截 距 长 度 平 均 密 度 /km'c 


垂直 TEC 平 均 误差 /TECu 


iE: a) GPS 地 基 观 测 数 据 水 平 穿 刺 信 息 ;b) 多 模 地 基 及 掩 星 观测 数据 水 平 穿刺 信息 ;c) 初始 场 垂 直 TEC 的 
平均 误差 ，d) GPS 地 基 观 测 数据 同化 结果 的 垂直 TEC 的 平均 误差 ，e) 多 模 地 基 及 掩 星 观测 数据 同化 结果 的 垂 
直 TEC 的 平均 误差 。 


图 9 日 均 水 平 穿刺 信息 及 垂直 TEC 的 平均 误差 分 布 


当 观 测 数据 增多 时 ， 观 测 数据 用 于 同化 过 程 中 需要 更 多 时 间 和 内 存 进 行 格 网 化 构建 
观测 算 子 ， 及 进行 方程 求解 。 许 多 工作 基于 超 算 使 用 公共 地 球 系 统 模型 (community earth 
system model, CESM) 或 DART (data assimilation eesearch testbed)， 对 电离 层 使 用 结构 格 
网 或 跨 尺 度 预测 模型 (the model for prediction across scales, MPAS) 非 结构 格 网 电离 层 模型 
进行 同化 分 析 ””; 还 有 许多 研究 人 员 通过 使 用 稀 疏 矩阵 及 稀疏 方程 求解 等 方式 降低 电离 
层 同 化 中 对 内 存 和 计算 量 的 极 大 要 求 ””" ; 为 得 到 三 维 电离 层 主要 变化 规律 或 者 异常 分 
布 特征 ， 基 于 非 均 匀 电 离 层 格 网 划分 构建 电离 层 同化 模型 也 可 以 在 保证 关注 区 域 精度 的 情 
况 下 降低 计算 需求 ””。 
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5 电离 层 应 用 及 展望 


通过 利用 理论 模型 /经 验 模型 与 多 源 电离 层 观测 数据 进行 融合 ， 获 取 精 细 三 维 电 离 层 
结构 及 变化 有 诸多 方面 的 应 用 与 研究 。 通 甬 过 对 流 层 、 热 层 、 等 离子 体 层 与 了 电离 层 耦 合 的 研 
究 呈 呈 拉 ， 可 以 揭示 电离 层 磁 暴 物 理 过 程 呈 汪 ， 暴 雨 时 期 闪电 在 电离 层 与 平流 层 的 作用 过 
程 ”， 厄 尔 尼 诺 -南方 涛 动 (the El Niño- Southern Oscillation, ENSO) 与 TEC 之 间 的 关 
系 ， 以 及 监测 电离 层 同 震 等 。 

电离 层 的 研究 依然 是 多 元 而 兴盛 的 领域 。 从 气象 数据 同化 在 物理 及 观测 数据 的 深度 耦 
合 ， 到 流体 力学 有 限 元 分 析 方法 可 实现 高 动态 流体 的 数值 模拟 收敛 ， 它 们 将 会 对 电离 层 数 据 
同化 有 很 大 的 指导 意义 ， 另 一 方面 ， 新 的 算法 与 新 的 算 力 的 加 入 ， 使 得 我 们 可 以 获得 更 加 精 
细 的 电离 层 时空 尺 度 特征 。 
5.1 ” 自 适 应 格 网 在 电离 层 同化 中 的 使 用 

电离 层 层 析 或 同化 需要 对 电离 层 空间 进行 格 网 化 。 可 依据 电离 层 中 电子 密度 的 高 度 分 
布 规律 ， 在 电离 层 密度 高 度 梯 度 变 化 较 大 区 域 使 用 较 细密 的 划分 ， 在 其 变化 较 小 的 区 域 使 用 
较 稀疏 的 划分 。 这 种 划分 方式 被 应 用 到 新 版 IRI Plas ^ US Jt Wu 等 人 加 对 电离 层 顶部 和 
等 离子 体 层 同 化 的 工作 中 ， 从 而 获得 电子 密度 垂直 分 布 上 的 细节 。MPAS 是 一 个 合作 项 目 ， 
旨 在 开发 用 于 气候 、 区 域 气候 和 天 气 研究 的 大 气 、 海 洋 以 及 其 他 地 球 系统 模拟 组 件 。 它 通 
过 非 结构 化 Voronoi 网 格 和 C 网 格 离散 化 ， 允 许 球体 的 准 均匀 离散 化 和 局 部 精细 化 ， 也 可 
以 预测 格 网 边缘 速度 的 正常 分 量 ， 特 别 适合 于 高 分 辨 率 、 中 尺度 大 气 和 海洋 模拟 一。 上 
述 两 种 方法 都 很 难 顾及 到 变化 的 电离 层 峰 值 高 度 在 电离 层 建 模 中 产生 的 影响 。 如 果 可 以 使 
有 限 元 方法 中 的 自 适 应 Mesh 算法 ， 对 于 电离 层 格 网 化 依照 电子 密度 的 空间 分 布 结 构 来 规 
划 ， 可 以 实现 对 电离 层 同 化 的 空间 分 辨 率 和 计算 时 间 的 优化 。 通 过 不 同时 刻 的 电子 密度 场 的 
细节 及 其 相关 性 来 控制 其 电离 层 相 对 偏 大 和 快速 变化 部 分 的 细节 ， 实 现 同 化 结果 的 三 维 相 
关 性 和 空间 整体 连续 性 。 
5.2 “机 器 学 习 在 电离 层 同化 中 的 使 用 

基于 多 种 电离 层 探测 手段 ， 人 们 通过 多 年 电离 层 观 测 收集 了 海量 多 源 电离 层 数据 。 深 度 


T di BARA rp Re. FRAMAR ARERR, KIRATA E A R A 
Q 识 ， 同 时 也 会 为 电离 层 异常 预报 与 危害 控制 提供 大 量 参考 信息 。 长 期 斜 向 电子 总 量 、GIM 
及 ionPrf 数据 等 通过 人 工 智能 算法 进行 电离 层 电子 密度 、 异 常 分 类 、 参 数 预测 或 者 GIM 分 


布 预 测 ， 都 已 陆续 开展 或 已 经 有 了 初步 研究 

当然 人 工 智 能 算法 不 仅仅 对 观测 数据 或 得 到 的 产品 (如 GIM 和 电子 密度 ) 等 进行 分 类 、 
拟 合 或 预测 ， 在 构建 三 维 电 离 层 模型 算法 中 使 用 深度 学 习 框 架 也 会 有 诸多 优点 : 在 电离 层 同 
化 算法 中 ， 通 过 深度 学 习 框 架 (如 tensorflow X: GraphBLAS"" ) 已 有 的 对 OpenCL (open 
computing language) #% CUDA (compute unified device architecture)” 的 集成 ， 将 使 得 
算法 计算 速度 得 到 极 大 的 提升 。 通 过 深度 学 习 框 架构 造 代价 函数 ， 然 后 利用 多 种 优化 策略 ， 
可 使 电离 层 同 化 算法 获得 极 大 改善 。 通 过 并 行 算 法 及 稀疏 矩阵 可 以 极 大 节省 物理 内 存 并 提 
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高 电离 层 数据 模拟 及 反 演 计算 效率 。 
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Current Status and Development 
of 3-D Ionospheric Research 
Based on Multi-source GNSS Observation Data 


FU Nai-feng?, CHANG Ying-li? 


(1. Shanghai Astronomical Observatory, Chinese Academy of Sciences, Shanghai 200030, China; 2. U- 
niversity of Chinese Academy of Sciences, Beijing 100049, China; 3. College of Information Technology, 
Shanghai Ocean University, Shanghai 201306, China) 


Abstract: In recent years, with the development of GNSS (Global Navigation Satellite 
System) technology, multi-source GNSS ionospheric observation data utilizes the spatial 
characteristics of ground-based and space-based GNSS ionospheric observations to form a 
more complete three-dimensional ionospheric observation system, promoted the progress of 
the ionosphere in theoretical models, empirical models, and assimilation models also pro- 
vided observations and methods for physical mechanisms of extreme weather, solar activity, 
and magnetic storm events by combining atmosphere and plasmasphere researches. The 


development status of ground-based and spaceborne GNSS ionospheric observation data, 


1 期 付 乃 锋 ， 等 ， 基 于 多 源 GNSS 观测 数据 的 三 维 电 离 层 研究 现状 及 发 展 81 


and the progress of spaceborne GNSS ionospheric observation data in the improvement of 
the empirical model of the topside ionosphere were displayed. The new developments and 
requirements in the process of an ionospheric assimilation model under the background of 
increasing multi-source GNSS ionospheric observation data are summarized, and it is point- 
ed out the following tips: (1) Due to the higher accuracy of the ionospheric empirical model, 
it can be used as the initial field of the assimilation model. (2) The theoretical ionospher- 
ic model can be used as the dynamic part of the Kalman filter. (3) It is possible to use 
multi-source data to reconstruct or predict the ionosphere with assimilation. And, it also 
pointed out that technologies such as adaptive meshing and machine learning can become 
new research directions of ionospheric assimilation algorithm to promote the research of 3D 


ionosphere. 


Key words: GNSS; ionosphere; empirical models; assimilation models 


